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Abstract 

The ground-state energy of a three-dimensional polaron gas in a magnetic field 
is investigated. An upper bound for the ground-state energy is derived within a 
variational approach which is based on a many-body generalization of the Lee-Low- 
Pines transformation. The basic contributing ingredients found are the ground-state 
energy and the static structure factor of the homogeneous electron gas in a magnetic 
field. Both these quantities are derived in the Hartree-Fock approximation. The 
resulting ground-state energy of the polaron gas is analyzed as a function of the 
electron density and of the magnetic field strength. 
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1 Introduction 



Increasing interest in the polaron gas problem is stimulated by intense ex- 
perimental investigations of polar materials (see, for review, Refs. [1,2]). The 
electron-phonon interaction is manifested, in particular, in experiments on the 
cyclotron resonance in bulk [5] and quasi-2D [6,7,8] systems, giving evidence 
for a resonant magneto-polaron effect [9] . Cyclotron-resonance measurements 
on semiconductor quantum wells with high electron density [10,11,12] reveal 
that the CR lines are shifted and split as compared to those in low-density 
quantum wells. In Ref. [13], these effects are interpreted within the framework 
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of the polaron concept, which thus remains vahd for high-electron density sys- 
tems. In this connection, the investigation of the ground-state properties of 
an interacting polaron gas is an important ingredient. 

In the present work, we analyze an interacting polaron gas at T = 0. However, 
the obtained results can also be applied to a polaron gas at nonzero tempera- 
ture if the thermal energy is small as compared to both the LO-phonon energy 
and the Fermi energy {ksT -C huji^o, ksT -C Ep). An upper bound for the 
ground-state energy of a polaron gas in an applied homogeneous magnetic 
field is derived using a many-body generalization [14] of the familiar Lee- Low- 
Pines transformation [15] for the polaron. The Hamiltonian of the polaron gas 
is given by 

^ = 2^ E / d'r'^M (7 V + ^ A(r) j ' *.(r) 

+ j dh'-^l{v)-^l{v')v{v-v')-^,,{v')^>,{v) 

+ E / d\^\{x)giiBBa^,{v). (1) 



where aj^ and Ok are the creation and annihilation operators for the longitudi- 
nal optical (LO) phonons with frequency c^loi 9 is the effective Lande factor, 
/^B is the Bohr magneton, m is the band mass, and ^'j^(r) and \E'cr(r) are the 
fermion creation and annihilation field operators for the charge carriers in po- 
sition r with spin projection cr = ±1/2 along the z axis which is chosen in the 
direction of the magnetic field B. The electron-phonon interaction strength is 
described by the amplitude 

(2) 

where a is the dimensionless electron-phonon coupling constant with the high- 
frequency and static dielectric constants and Eq, and v{r) is the Coulomb 
interaction between the electrons, with the Fourier transform 

v{k) = (3) 

where V is the crystal volume. 

In Rcf. [14], the ground-state energy was studied for a polaron gas in the ab- 
sence of a magnetic field by introducing the following unitary phonon transla- 
tion operator, which is the many-particle generalization of the Lee-Low-Pines 
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transformation [15]: 



U — exp 



2:/.v*i(.)(i:(w--/;4e-))*„w 



(1) 



By applying exactly the same unitary transformation for the polaron gas 
in a the presence of a magnetic field and calculating the expectation value 
{^GT \U~^HU\ ^gt) in a trial ground state \^gt)= \vac) ■ |\l/ei) which is the 
product of the phonon vacuum and an electron state, we found an upper bound 
to the ground-state energy 

k 

where Eb is the energy of a gas of non-interacting electrons in an applied 
magnetic field B, and >S'(k) is the static structure factor of a homogeneous 
electron gas in a magnetic field. Minimizing (5) with respect to we obtained 
the same formal expression for the optimal values of /k as in [14] : 



- 7 : — Wk^- (6) 



This results in the following expression for the variational ground-state energy 
per particle: 

Eo-Eb- W^k) [1-S (k)] - E ^^'^[^^^l ■ (7) 



The ground-state energy (7) is structurally similar to that obtained in [14] in 
the absence of a magnetic field. The basic ingredients Eb and >S'(k) are different 
though, because of the magnetic field. The static and dynamic structure factors 
can be calculated using various approximations (see, e. g., Ref. [16] for a 
confined 3D electron gas and Refs. [17,18,19] for a 2D electron gas). 

In the present paper, we calculate Eb and S'(k) in the Hartree-Fock approx- 
imation, which is variational in nature. The contribution Eb is treated in 
Section 2.1. The static structure factor S{k) will be discussed in Section 2.2. 
The upper bound to the ground-state energy per particle, taking both contri- 
butions into account, is discussed in Section 3. For the numerical calculations 
and illustrations the degenerate polar semiconductor GaAs is used (Table 1). 
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Table 1. Pcirameters and material constants of GaAs which are 

used for the calculations. 



m/me 


tuji^o (meV) 






a 


9 


0.0653 [20] 


36.3 [21] 


10.89 [21] 


13.18 [21] 


0.068 [22] 


-0.44 [23] 



2 Ingredients contributing to the ground-state energy 



2.1 Energy of non-interacting electrons in a magnetic field 



The energy of non-interacting electrons in a magnetic field is calculated start- 
ing from the energy for a single particle in a magnetic field and applying the 
Pauli principle. Expressing that the number of particles N equals the number 
of occupied states we arrive at the equation 



(8) 



where v is the degeneration factor of each Landau level, kz represents the wave 
vector component along the magnetic field, uOc is the cyclotron frequency and 
Ep stands for the Fermi energy. The function (x) is defined as 



1 if X is true 
if a; is false 



(9) 



By summing all one-electron energies obeying the Pauli principle one obtains 



Ep 



mwr 



27r^^no 



»imax(a-) 

E E 



^2 ,^3 

6m 



(10) 



«n,a = Y ^ (-^^ -hu;c(n + ^^ - giJiBBa^ , (11) 

where uq = N/V is the electron density and K.n,a is the maximal wave vector 
in z-direction of an electron with the spin projection a in the n-th Landau 
level. The integer 

'Ep — qubBg r 
hujf. 2 . 

denotes the number of fully occupied Landau levels with the spin projection 
(7. The energy contributions clearly depend on the electron density and on the 



(12) 
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magnitude of the applied magnetic field via the Fermi energy and the number 
of occupied Landau levels. 



2.2 Structure factor 



Within the Hartree-Fock approximation, the expression for the static structure 
factor is calculated based on the second quantization form 

era' 

X {^,1 (r)^^,(r')^,,(r')^,(r) . (13) 



The calculations were performed with periodic boundary conditions on a cube 
with volume V, and the limit K — > oo is taken keeping the density tiq con- 
stant. Writing out the field operators in the energy representation and ap- 
plying Wick's theorem, the structure factor is expressed in terms of integrals 
over wave vectors and summations over eigenstates of electrons in a magnetic 
field. Using the orthogonahty relations for the eigenstates one is left with an 
exchange term which involves integrals of the form [24] 

oo 

J dx e-'''Hm(x+a)Hr,(x+(3) = 2'*V7rm!/3"-"*L;^-"*(-2Q;/3) for m < n (14) 

— oo 

where Hn{x) is the Hermite polynomial of the n-th order and where Lj{x) is 
an Laguerre polynomial. 

The resulting expression for the static structure factor of an electron gas in a 
magnetic field is 

,,,, , 1 'l^i -..ax(<x)n M , /;|^2 X n>-n< 



X 



rn>-n< / ^B^J 



C{n,n',a,k^). (15) 



^{n, n , (7, kz) — 2Kjj^^(^0 (l^^l ^ l^n<,a l^n>,(j) 

~l~ {^n^,a ~l~ ^n^,a |^z|)Q {^n^,a ^n^,a ^ |^z| ^ ^ny,a ~l~ ^n<c,cr^ 

(16) 

where Ib = \[^ is the magnetic length, k_[_ denotes the component of the 
wave vector perpendicular to the applied magnetic field, n> = max(n, n') 
and n< = min(n, n'). In order to avoid numerical underfiow or overfiow in 
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the recurrence relations, it is important to consider renormalized associated 
Laguerre polynomials 

np / X _ J p\{n-p)\ 



It is clear from (15) that the dependence of 5'(k) on the wave vector component 
perpendicular to the magnetic field and on the wave vector component along 
the magnetic field can be treated separately, and that the introduction of a 
magnetic field results in an anisotropy of the structure factor. 

In the present calculations, the electron density is described through the pa- 
rameter Ts = ro/a%, where tq = (47mo/3) is the Wigner-Seitz radius, and 
a*^ = Y?Zool (e^rn) is the effective Bohr radius. 

In Fig. 1, the static structure factor S {k±, kz) is plotted for the value = 1, 
corresponding to the electron density no ~ 3.75 x lO^^cm"^ in GaAs. The 
dependence of the structure factor on k±_ is shown in Fig. 1(a), which reveals 
that the structure factor S {k^,kz)\^^^Q is a smooth function with limiting 
values jS" (0, 0) = and (oo, 0) = 1. The curvature however increases with 
increasing magnetic field, with a smaller slope for small /c, and an increasing 
value of the wave vector before the high wave vector limit is reached. 




k^kp ^J^F 



Fig. 1. Hartree-Fock static structure factor S {k±,kz) of an electron g function 
of the wavevector component k± perpendicular to the applied magnetic field (a) and 
as a function of the wavevector component kz along the applied magnetic field (b). 
In both cases the other wavevector component is equal to zero. In the panel "b", 
the downward pointing arrows indicate the two kinks corresponding to the magnetic 
field of 4 Tesla (nmax = 2), and the horizontal pointing arrow indicates the kink 
corresponding to 7 Tesla (nmax = !)■ 

Fig. 1(b) reveals that the structure factor as a function of the wave vector 
component along the magnetic field is far less smooth. Relatively pronounced 
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kinks appear in S (0, kz) where the slope suddenly changes, and the number 
of kinks for S (k) < 1 exactly equals the number of fully occupied Landau 
levels. The static structure factor S{k) thus allows to count the number of 
occupied Landau levels of a polaron gas with given density and magnetic 
field strength. Furthermore, S {k±, kz) more rapidly tends to unity in the z- 
direction, parallel to B, than in the xy-plane. Furthermore, this anisotropy 
enhances with increasing B. As a result, S {k±,0) general trend, a 

decreasing function of B, whereas S (0, k^) increases with increasing B. The 
physical reason for this anisotropy is the following. A magnetic field applied 
along the 2;-axis hinders the translation of the electrons in the xy-plane anal- 
ogously to a parabolic confinement potential with the frequency cUc and with 
confinement length Ib = \f^- The magnetic length decreases with increasing 
magnetic field strength, and hence the region of k±, in which S {k±,kz) sub- 
stantially varies, increases because it is of the order of as seen from the 

factor exp ( — ) in Eq. (15). However, the anisotropy is non- monotonous: 

it is accompanied by oscillations due to the de Haas - van Alphen effect [25]. 



3 Ground state energy: results and discussion 



Converting the sums in (7) to integrals over cylinder coordinates, introducing 
dimensionless variables, reducing the singularity in the origin and using two- 
dimensional open Romberg integration, we calculated the variational ground- 
state energy per particle numerically. 
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Fig. 2. The ground-state energy per particle (a) and the Fermi energy (b) as a 
function of the magnetic field for a polaron gas in GaAs with = 1. The arrows and 
numbers indicate the changes of the number of filled Landau levels nmax- The dashed 
curve shows the ground-state energy calculated without the polaron contribution. 
Notice the different scale in the abscises of both panels. 

In Fig. 2(a) the ground-state energy per particle is plotted as a function of 
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the magnetic field. Both the many-polaron ground-state energy and the Fermi 
energy are oscillating functions of the magnetic field. At well-defined magnetic 
fields infiection points occur in Eq (B) . This phenomenon can be understood 
if we look at the Fermi energy: the inflection points of the ground-state energy 
coincide with the local maxima of the Fermi energy, and therefore those inflec- 
tion points are a direct consequence of the de Haas - van Alphen effect. For 
comparison, we also have plotted in Fig. 2(a) the ground-state energy calcu- 
lated for q; = (the dashed curve), i. e., without the polaron contribution. As 
seen from the figure, the electron-phonon interaction substantially shifts the 
ground-state energy to lower values with respect to that calculated for a — 0. 
This shift is practically independent on the magnetic field strength. 

The described dependence of the energy on the magnetic field is also seen in 
Fig. 3(a), which shows the energy as a function of the parameter r^. From 
this figure, it is clear that by applying a magnetic field, the minimum of the 
ground-state energy occurs at a finite value of r^. 
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Fig. 3. The ground-state energy per particle with and without applied magnetic field 
(a) and the Fermi energy for 5 = 4 T (b) for a polaron gas in GaAs as a function of 
the parameter Vg- The arrows and numbers indicate the changes of rimax- The dotted 
curve shows the ground-state energy calculated without the polaron contribution. 

Similarly as in Fig. 2(a), kinks appear due to the change of the number of 
fully occupied Landau levels, as is more clearly illustrated in Fig. 3(b) where 
the Fermi energy is shown as a function of r^. For all values of the magnetic 
field, the magnitude of the difference between E'^ and falls down with 

increasing density, so that the influence of the magnetic field on the ground- 
state energy of a polaron gas is larger for lower densities. This behavior of E'^ 
as a function of the density is explained by the fact that, when the density 
is high enough, a large number of Landau levels are flUed. In this case, the 
magnetic field should be considered as relatively weak with respect to the 
characteristic frequency of the electron gas Ep / h {lUc -C Ep / h) . 

When comparing in Fig. 3 the density dependence of the ground-state energy 
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for a polaron gas with that calculated without the electron-phonon interaction, 
we see that the magnitude of the polaron contribution decreases with decreas- 
ing Ts due to screening. The same trend was indicated also for a polaron gas 
in the absence of a magnetic field in Ref. [14]. 



4 Conclusion 

We have extended the variational approach to the many-polaron problem [14] 
to a system of interacting polarons in an external homogeneous magnetic field. 
The rigorous upper bound for the ground-state energy of an interacting po- 
laron gas in a magnetic field has been obtained in terms of the static structure 
factor. Due to the presence of a magnetic field, the static structure factor as 
a function of the electron density reveals features related to the fiUing of the 
Landau levels, when the Fermi energy crosses a Landau level. The other result 
of an applied magnetic field is the field-induced anisotropy of the structure 
factor S (k) , which should be observable experimentally. 

The apphcation of a magnetic field has pronounced effects on the behavior 
of the ground-state energy as a function of the electron density. The most 
prominent feature, which is also observed in the ground-state energy as a 
function of the magnetic field, is the occurrence of kinks in the ground-state 
energy, which occur as a consequence of the de Haas - van Alphen effect as soon 
as a Landau level becomes fully occupied. Although the ground-state energy is 
not directly observable, these phenomena should infiuence experimental data, 
e.g., the optical absorption of a polaron gas in a magnetic field, which are 
currently being studied. 
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